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Organic Electrochemical transistors (OECTs) present unique features for their strategic combination with 
biomedical interfaces, simple and low voltage operation regime and sensing ability in aqueous environment, 
but they still lack selectivity, so that a significant effort in research is devoted to overcome this limitation. 
Here, we focus on the diffusion properties of molecular species in the electrolyte, which opportunely 
analyzed, modeled and compared to experimental data, serve as a simple and direct key factor in the 
recognition of species during OECT sensing. Specifically, we model the transient behavior of an OECT 
considering the effect of diffusion of the target species in the electrolyte. In doing so, we develop a general 
method that can be used to differentiate and distinguish different molecules from a complex mixture, on the 
basis of their diffusivity and thus mass. More importantly, the model can be realistically used to determine 
the physical characteristics of the transported species in a solution from a simple fitting procedure. On the 
basis of the obtained results, we discuss the contribution that our study could bring to OECT architecture to 
realize a new generation of devices with improved sensitivity, selectivity and reliability. 

Since the discovery of conducting polymers in the mid-1970s 1 , the development of innovative devices based 
on their structural and conducting properties has grown interest for application in organic electronics and 
bioelectronics. Organic Field effect Transistors based on different organic structure have been realized, with 
promising application in printed and flexible electronics 2 . Among the different approaches to face organic 
electronic challenges, much attention has been dedicated to devices which exploit the interactions of organic 
polymer with ionic electrolytes: Organic Electrochemical transistors (OECTs) present unique features for their 
strategic combination with biomedical interfaces 3 , simple and low voltage operation regime 4 and their sensing 
ability in aqueous environment 5 . OECT bases its working principles on the de-doping action of a liquid or gel 
electrolyte over the conducting channel of a polymeric film. The interaction, between ionic species in the 
electrolyte and the active area of the organic conductive polymer, controls the modulation process in the 
transistor, used for either electronic or sensing applications 6 ' 7 . One of the most used p-type organic semi- 
conductor is Poly(3,4-ethylenedioxythiophene) doped with poly(styrenesulfonate) (PEDOT:PSS), which pre- 
sents high stability in a wide pH range, high conductivity, low cost and lithographic process integration 8 . OECT 
devices based on PEDOT:PSS have been employed in different technological fields such as chemical and bio- 
logical sensing 9 , control of cell adhesion 10 and viability 11 , sport and healthcare monitoring through integration on 
natural cotton fiber 12 . Among the different applications, one of the crucial challenges to realize effective iden- 
tification of the ionic species in-sensing and bio-sensing is to develop and improve electivity properties. The 
ability to reveal with high sensitivity different species in electrolyte solutions has been demonstrated 713 , never- 
theless recognizing the different species revealed is the actual challenge to be faced, and this requires different 
strategies 1415 . Among the possible approaches, one of the most simple and direct is to evaluate the different 
diffusion properties of the ionic species in the electrolyte. In considering the physical characteristics which 
influence the diffusion properties of the ionic species it is possible, using a correct model for the OECT working 
processes, to distinguish different ionic species. Previously, an important model for charge transport in organic 
electrochemical transistor has been proposed by Bernards and al. 16 , which takes care of a complex interplay 
between ionic and electronic motion. In particular, a separation between ionic transport in the electrolyte and 
hole transport in the PEDOT:PSS channel is considered. The result is that the hole electronic charge is modeled as 
a semiconductor, while the ionic charge is considered on the direct effect of a drift voltage. Using an ionic circuit 
equivalent toideal polarizable electrodes, this model reproduces the steady stateand transient characteristics of 
PEDOT:PSS OECTs, but it disregards the effect of diffusion of the ionized species in the electrolyte. 

Here, we revised the model of Bernards to consider the effectof diffusion. Using the formalism of the 
Nernst-Planck equations 17 " 19 (Methods) we described the motion of the charged species in the electrolytein terms 
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of the diffusivity and charge of those species (Figure 1). This motion, 
in turn, results crucial to determine the effective behavior of drain 
source current in dependence to the physical properties of ionic 
species. 

In comparing experimental results of a real OECT device with the 
predictions of our theoretical model, we recognized and identified 
four different species (that are, Na, Mg, K, Zn) in a solution on the 
basis of their different mass and charge (and thus diffusivity). 
Moreover, from a simple fitting procedure, we estimated the numer- 
ical values for the diffusion coefficients of the released species with 
good accuracy. 

The model effectively predicts the different kinetic behavior of the 
source-drain current for the different species, allowing a recognition 
and identifications of the ionic salts, improving the selectivity of the 
device. The results pave the way to an innovative kinetic analysis of 
the OECT modulations, useful for selective recognition of ionic spe- 
cies, with many sensing and bio-sensing applications. 

Results 

Effect of diffusion. On varying the diffusion coefficient over a 
significant range, we observed a very high sensitivity of the tran- 
sient behavior of the OCET on the diffusivity of the ionized species 
in the electrolyte. Fig. 2a reproduces the simulated transient current 
I(t) as a function of time for the effective diffusion coefficient varying 
from 10" 1 to 10 times D G , where D G = 10" 8 m 2 /s is an intermediate 
value of diffusivity used as a reference. In the considered case, z = 1 is 
hold fixed, while the remaining variables utilized in the simulations 
were adjusted as follows: (p = 1 V, x = 1 cm, Ex = — d(p/dx = 100 V/m, 
T = 298 K, jll = 5 10 26 m 2 /(V s), p G = 5 10" 4 m" 3 eV"\ q = 10" 19 C. 

For fixed values of D, the current increases for increasing time; 
however, the exact time behavior of the system would significantly 
depend on the magnitude of D. For low values of D, that is, for D/D G 
< 1.5, the sensor behaves as if the output conforms to the standard 
curve of a first- order response, where the rise time and the time 
constant, that is here defined as the time which is required by the 
sensor to reach 90% of the steady-state or maximum level upon 
application of the external voltage, diminish with D. For values of 
D/D 0 larger than 1.5, the time-dependent characteristic of the system 
resembles that of a slightly undamped second order system, here, the 
response is even faster with the typical overshoot of similar dynamic 
systems. Fig. 2c reports the simulated current as a function of y = D/ 
D Q for defined values of time, the derivative dl/dy of this diagram is 
reported separately in Fig. 2e; they both provide information about 
the sensitivity of the current to the diffusion coefficient. Collectively, 
these results indicate that the effect of diffusion is moderate at the 
early time of the process, while it is significant for larger values of 
time. Moreover, the sensitivity of the sensor would depend on D. For 
vanishingly small values of diffusion, that is, for y —> 0, the intensity 
of current will vary significantly for small variations of D. At any time 
of the simulated process of measurement, the current is proportional 
to D, and this is easily explained considering that augmenting D is 



equivalent to using smaller molecules, that can be more easily trans- 
ported from the electrolyte to the PEDOT conductive polymer. 

Effect of charge. In this section, the effect of charge on the overall 
response of the system is analyzed. In doing so, the diffusion 
coefficient is maintained constant as D = Do, while the charge 
number z is varied from 0.1 to 10. The simulated transient current 
I(t) is reproduced in Fig. 2b as a function of time for these prescribed 
values of z. For a fixed z, the current smoothly transitions from zero 
to a steady state value that increases for increasing z, and the sensor 
behaves like a first order system. Differently from the case of 
diffusion, the behavior of the system is similar for any charge 
number, and thus z may be considered like a simple scale factor. 
This is even more evident in Fig. 2d, where I is reported as an 
explicit function of z at different times. When t is small, that is, 
approximately for t < 8 s, the current I is constant over z. For 
large values of time, that is, t > 8 s, I is quasi linear in z, and thus 
the effects of z on I are easily predicted. However, these effects are 
small, in the sense that the variations of I are limited on changing z 
(Fig. 2d). The analysis of sensitivity to z (Fig. 2f) confirms this 
understanding. Overall, on comparing Fig. 2e and Fig. 2f it appears 
that the contribution of diffusion on the absolute magnitude of ionic 
current dominates over the contribution of charge number. Or, that 
is equivalent, the system responds more promptly to variations in D, 
compared to variations in the charge number z. 

Combined effect of diffusion and charge on the current settling 
and rise time. The model here introduced makes it possible to 
predict the current response of an OECT as a function of diffusion 
and charge. The concerted action of these two parameters and their 
influence on the time constant! and rising time of the device is 
analyzed and reported in Figure 3. In Fig. 3a, the current settling 
time is shown in the form of a contour plot as a function of D/D 0 
and z/z G (D G = 10" 8 m 2 /s, z G = 1). The time to the steady state 
behavior of the device is large for small values of the diffusion 
coefficient, and gradually decreases for increasing D. Notice, from 
the diagram, that D dominates over z. For vanishingly small values of 
D to nearly D = 4D G , the influence of z can be neglected, and the time 
constant does not vary on changing z. This scenario changes for 
larger values of D. On augmenting D or, equivalently, on 
considering smaller molecules, the effect of the charge cannot be 
disregarded. Low values of t are thus obtained in the region of the 
diagram where D and z are simultaneously large. For the particular 
set of values considered here, the time constant of the systems varies 
from —60 s, corresponding to diffusion coefficients smaller or equal 
than about D G , to less than 13 s, that instead is obtained for D and z 
larger than, respectively, 8 Do and 3 z G . 

An important parameter for characterizing the time response of 
the system, is the rising time k, that is here defined as the time 
derivative of the ionic current calculated at the early time of the 
process (that is, t = 0). It gives a measure of how rapidly the current 
increases. In Fig. 3b, k is reported as a function of D and z. Similarly 




Figure 1 | Schematic representation of the ionized species within the electrolyte that, upon the application of an external voltage, are driven towards the 
conductive polymeric PEDOT:PSS film. The motion of the ions depend on pure drift(that is, the constant electric field) and pure diffusion. 
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Figure 2 | The simulated transient current I(t) as a function of time for the effective diffusion coefficient varying from 10 _1 to 10 times D Q , where 
D Q = 10~ 8 m 2 /s is an intermediate value of diffusivity used as a reference (a). The simulated transient current I(t) as a function of time for the charge 
number z varying from 0.1 to 10 (b). The simulated current as a function of y = D/D 0 for defined values of time (c).The derivative dl/dy, that is, the 
sensitivity analysis of I to D (e). I is reported as an explicit function of z at different times (d). The analysis of sensitivity of I to z (f). 



to the other described case (effect on time constant), the rise time of 
the current is dominated by diffusion for small values of D. The 
influence of z emerges for D larger than about 3 D G . Interestingly, 
the values of k range from —0.04 to —0.4 for D varying from 0.5 to 
about 9 times D, and z from 1 to 4. That is to say that a tenfold 
increment in k is obtained for a combined increment of D and z of 
nearly two orders of magnitude. 



Figures 3 a and 3b are very similar to the concept of design maps for 
OECT sensors. Somewhat similar in concept to the periodic table of 
chemical elements, but more akin to phase diagrams in their graphic 
presentation, these design maps may provide guidance for the 
development of optimized OECT-biosensors and are for the time 
being limited to the combination of diffusion and charge. These 
maps and their more sophisticated evolution that will be developed 
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over time will allow for detection regimens to be personalized to 
individual cases, resulting in greater efficacy and sensitivity of the 
measurement and analytes detection. 

Determining the diffusion coefficient of unknown species in a 
solution. The described model is predictive in nature as 
demonstrated in the above described results. More importantly, it 
may be realistically used to determine the physical characteristics of 
the transported species in a solution. This can be realized on 
inverting the model on experimental data. To demonstrate this, 
four different ions were considered in a real OECT device as 
described in the Materials and Methods, that are, Magnesium 
(Mg ++ ), Zinc (Zn ++ ), Potassium (K + ), Sodium (Na + ). These 



elements feature different diffusion coefficients, ranging from D = 
0.706 10" 9 m 2 /s and D = 0.703 10" 9 m 2 /s for Magnesium and Zinc, 
to D = 1.334 10" 9 m 2 /s and D = 1.960 10" 9 m 2 /s, for Sodium and 
Potassium 20,21 . Moreover, the charge number of those diffusion 
species is also different, being z = 1 for sodium and potassium, 
and z = 2 for magnesium and zinc. Therefore, one would expect a 
different time dependent behavior of the device on changing the 
probed species. This is indeed the case for the measured ionic 
current. In Figure 4 the experimental current sensed by the device 
is reported as a function of time for the considered species in 
solution. The transient behavior for Zn and Mg is similar in shape 
and this is easily accepted considering that Zinc and Magnesium 
have physical constants that are comparable. In the same diagram, 
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Figure 4 | Experimental current sensed by the OECT device reported as a function of time for four different species in solution: Mg, Zn, Na, K. 

The experimental data were fitted using the described diffusion model; from this fitting procedure, the numerical values for the diffusion coefficients of 
the released species were estimated with good accuracy. 
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Figure 5 | On patterning the PEDOT polymeric film at the nano scale, an increasing degree of complexity is introduced into the geometry of the 
devices. This extra not-continuous scale would provide the device with additional functions, including the capability of manipulating extremely small 
volumes of a solution. For this configuration, the motion of charged species in the electrolyte is determined from the interplay of convection, diffusion, 
and evaporation. 



differently from the case of Zn and Mg, and with a lower intensity, the 
ionic current of Potassium (in light red) and Sodium (in light blue) 
are reported. Notice that Mg, Zn, K and Na are ordered for decreas- 
ing intensity in the diagram in Figure 4, and this is consistent with the 
heretofore derived results. These sets of data were therefore fitted 
using our described diffusion model, from this fitting procedure, the 
numerical values for the diffusion coefficients of the released species 
were estimated as D Mg = 0.685 10" 9 m 2 /s, D Zn = 0.763 10" 9 m 2 /s, 
D K = 2.10 10" 9 m 2 /s, D Na = 1.213 10" 9 m 2 /s, that are similar to the 
true values of diffusivity, with errors, with respect to these, of e Mg = 
0.021, e Zn = 0,085, e K = 0.07, e Na = 0.09, that are small. This would 
demonstrate the ability of the model to infer the characteristics of a 
specific specie in a solution from a simple ion current measurement. 
Moreover, it demonstrates that the method may be utilized to 
differentiate different species in a complex mixture. 

Discussion 

The described results demonstrate that, in considering OECT sensors, 
the effect of diffusion and charge of the probed species in the electro- 
lyte cannot be neglected. In similar devices, the measurement of the 
ionic current, that is, the response of the system, results from a 
sequence of tightly interweaving processes, that are, (i) the transport 
of the species from the electrolyte into the conductive polymer, (ii) the 
injection of a cation into the organic film and the removal of a hole at 
the source electrode, (iii) the transport of holes between source and 
drain. Our results indicate that the onset of this chain of events, that is, 
(i) the transport of the species in the electrolyte, determines the final 
response of the system. The mechanisms of ion transportation in the 
electrolyte would depend, in turn, on the externally applied voltage, 
and on the electrical mobility of those ions in the solution, that is 
described by the well known Einstein-Smoluchowski relation 22 m = 
q/KT D. For a fixed temperature, the combined effect of the charge q 
and the diffusion D would therefore have an effective impact on the 



output of the sensor. Considering that the Brownian coefficient of 
diffusion is inversely related to the hydrodynamic diameter of a dif- 
fusing molecule, different species in a solution, with different sizes and 
thus mobilities, would behave differently upon the application of the 
same stimulus. The described mathematical method delivers the abil- 
ity to predict and anticipate those differences. This understanding may 
have some very practical applications in bio -sensing and molecular 
diagnostic. In the results, we demonstrated how this scheme can be 
utilized for detecting and recognizing specific species in a solution. 
More importantly than this, we demonstrated how, using a simple 
fitting procedure, one may retrieve the diffusion coefficients (and thus 
the atomic radii/masses) of unknown species in a solution, that is to 
say, the method is capable to perform the label free analysis of the 
analytical content of a solution on a quantitative basis. Its more soph- 
isticated evolution that will be developed over time will allow the 
detection of single species from a complex mixture. 

However, this concept has a certain relevance beyond these spe- 
cific results. Consider, for ease of visualization, the scheme in 
Figure 5. The theoretical model indicates that the displacement of 
the active ions in the electrolyte would ultimately determine the 
response of the system. For the present configuration, this displace- 
ment is driven by the externally applied electric fieldand by pure 
diffusion. It may be convenient, for specific purposes, introducing 
a supplementary parameter in the system i.e. the velocity vector field 
in the solution domain, an additional convective term. This can be 
done utilizing two different approaches. The first, utilizes active 
micro- and nano-fluidics, where active components as micro-pumps 
or micro valves manipulate extremely small volumes of fluids in 
constrained geometries at the micro and nano scale. This approach 
shall not be discussed here. The second (passive) approach, is 
based on an accurate texturing of the interface between the electro- 
lyte and the organic polymeric film. On patterning the PEDOT poly- 
meric film at the nano scale, an increasing degree of complexity is 
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introduced into the geometry of the devices. This extra not-continu- 
ous scale would provide the device with additional functions, includ- 
ing the capability of manipulating extremely small volumes of a 
solution. Depending on texture and chemical modification, a milli- 
metric drop would maintain a spherical shape if positioned upon 
those surfaces. This concept is at the base of super-hydrophobic 
surfaces, that are bio-mimetic, structured micro and nano-materials 
which reveal superior properties compared to conventional flat sur- 
faces, including extremely low friction coefficients 23 ' 24 . 

An aspect in regards to which a drop is different from the same 
volume of liquid constrained in a rectangular geometry, is that the 
drop develops convective flows within its volume due to the cele- 
brated Marangoni effect 25 . These flows of masses depend on the size 
of the dropand on its curvature, that is, on the angle of contact on the 
surface. Consequently, they may be predicted/calculated by control- 
ling the geometry of the substrate. Nanotechnologies provide a tool 
to fabricate the substrate on the basis of a specific design to obtain a 
convective flux with the desired characteristics. These flows repres- 
ent the additional convective term that may couple to diffusion to 
yield a complex field of motion that would ultimately modify the 
transportation of molecules from the electrolyte to the electrode, and 
thus the output of the measurement system. Where proven inad- 
equate, pure analytical models can be paired to computational fluid 
dynamics (cfd) methods and computer simulations to resolve the 
field of motion in complex geometries. 

Therefore, this scheme can be utilized to realize a system where 
diffusion and convection are separately controlled and, from the 
balance of these, different molecules in a mixture can be fragmen- 
tized, separated or selected. A correct design of the pattern can render 
this selection highly precise and reproducible. In doing so, a new 
generation of devices will be obtained with high sensitivity, selectivity 
and reliability. This concept has been exploited for some exploratory 
experiments reported in 26 , where a PEDOT conductive polymer and 
super-hydrophobic surfaces are, for the first time, merged together. 
For guaranteeing an efficient analysis, the fabrication of such devices 
should be always accompanied by suitable mathematical models, like 
that introduced in this or other similar works 27 ' 28 . The use of PEDOT 
has heretofore been limited to direct applications, if combined with 
nano -geometry, it could instead represent the infrastructureof a new 
devices which feature new, unprecedented functions. 

Methods 

Following 16 , the behavior of the OECT device is here analyzed on dividing the system 
into an electronic and an ionic circuit. The electronic circuit consists of ap-type 
organic semiconductor film that transports holes between source and drain electrodes 
and is described by Ohm's law. Differently from this, the ionic circuit accounts for 
transport of ionic charge in the electrolyte. The interaction of these two components, 
namely the transport of ions from the electrolyte into the organic semiconductor film, 
is fundamental to the behavior of OECTs. In a separated experimental section, we 
shall describe the experimental set up utilized to derive the data that confirmed the 
feasibility of the model. 

Transporting charged species into the channel. In what follows, we will describe the 
ionic transport of electroactive species upon the application of an external voltage. In 
doing this, consider the scheme in the Supporting Information Fig. SI. Assume that, at 
the initial time t = 0, n charges are uniformly distributed in a rectangular channel. 
Upon the application of an external voltage T, those charged particles will be 
transported across the channel. The driving forces responsible for this motion are pure 
diffusion and the electromagnetic force field. For this configuration, the distribution of 
charged particles can be described by the celebrated Nernst-Planck equation 17-19 



d( v«A+|!^V(p 



(1) 



that, for mono dimensional geometries described by the sole spatial variable x, takes the 
form 



f d\b zF , dm 



(2) 



in (2), J x is the molar flux, D is the Brownian or molecular diffusion coefficient, \|/ is the 
concentration, defined as the number of moles of charged particles per unit area, z is 
the charge number, F = 9.64 10 4 C mol" 1 is the Faraday constant, R = 8.314 J 



K" 1 mol" 1 the universal gas constant, T, the absolute temperature of the system; (p, 
instead, is the electric potential. We make now the assumption that cp is linear along the 
channel, being dty/dx = T = — E x , where E x is longitudinal component of the electric 
field, and thus (2) may be rephrased as 



-D 



zF 
RT 



using in (3) the continuity condition that guarantees the conservation of mass 
one would obtain the following governing equation 



or, that is equivalent 



S ± =D 

dt 



zF #' 

rt x a x 



dx 2 



(3) 



(4) 



(5) 



(6) 



and, on solving Eq.(6), the distribution ^(x^of the charged particles as a function of 
the spatial and time variables, would be determined. Eq.(6) has the form of a one 
dimension convection-diffusion solute transport equation, where the convective term v 
= D E x zF/RT is maintained constant for the extension of the channel. For the present 
configuration, the initial and boundary conditions are imposed as follows 



\J/(x,0) = \J/ o 0<x<l 
i^(oo,0=0 t>0 

Equation (6), along with the initial and boundary conditions (7), has an analytical 
solution that reads as 29 



(7) 
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2 J \2VDtJ 
where 1 is the channel length, and erf is the error function, defined as 

2 [* 2 



t Jo 



(9) 



the concentration in Equation (8) indicates how, for a fixed longitudinal distance x, the 
number of ions per unit area varies as a function of time. The total amount of charge 
that is transported into the OECT, that is, the PEDOT channel, may be accordingly 
derived as 



Q(t)=N A 



qS f 

Jo 



v\J/(1,t)-D 



# I 

dx \x=i 



dx 



(10) 



and Q(t) is again a function of time. In (10), N A = 6.02 10 23 is the Avogrado's number, 
q is the elementary charge, S is the cross sectional area of the Electrolyte. Q(t) would 
explicitly depend on D. 

Transient behavior of the drain current. The transient behavior of OECTs was 
solved in 16 .The problem was considered as the superposition of effects, that are, 
injection of a cation from the electrolyte into the organic film and removal of a hole at 
the source electrode. On neglecting the spatial variation of voltage and hole density, 
an expression for transient source-drain current I(t) was derived, being 



I(t) = G 1 



Q(t) 

qp„V 



V d -f 



gQW 

dt 



(11) 



where G = q up D V is the conductance of the organic semiconductor film, [i is the hole 
mobility, p D is the initial hole density, V is the volume of the semiconductor material, 
V d is the constant drain voltage, f is a corrective constant that accounts for the non 
uniform distribution of the hole density. Using in (1 1) Eq.(10), the transient current 
can be accordingly determined. Similarly to the total charge Q(t), I(t) depends on the 
molecular diffusion coefficient D. This is shown in Fig. S2, where I(t) is reported as a 
function of time for different values of diffusivity. 

Measuring the electric characteristic of ionized species in solution. Figure S3, in the 
Supporting Information, shows a sketch of the OECT realized on cotton fiber and of 
its electrical circuit. Figure S3 reports also a schematic representation of the 
mechanism underling the OECT sensing process. The saline species (NaCl, KC1, 
MgCl and ZnCl) disgregate in CI" and M + in water. In presence of a gate voltage, CI" 
ions reacts at the Ag electrode by an oxi-reduction process, while M + ions are forced 
to PEDOT:PSS channel. Hence a faradaic current, flowing in the source-gate circuit, is 
generated 5 . Under this OECT faradaic regime of operation, the potential drop at the 
electrolyte/ gate interface decreases and the effective gate voltage (V gse ff) increases, the 
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latter forcing H + cations to move toward the polymer surface and dedope 
PEDOT:PSS. 

Real time sensing of OECT were acquired by measuring Ids vs. time, fixing Vds at 
— 0.2 V and switching on Vgs from 0 V to 0.8 V at a fixed time, maintaining 0.8 V of 
gate potential for 60 s. The measurements have been performed for the different salts 
(NaCl, KC1, MgCl and ZnCl) with different masses and diffusion coefficient, with a 
concentration of 10~ 3 M. The application of Vds induces a drift of the holes along the 
PEDOT:PSS channel, generating a drain-source current (Ids). Upon application of a 
positive gate voltage (Vgs), cations (M + ) from the electrolyte enter the PEDOT:PSS 
channel causing its de-doping according to the following equation 4 : 



PEDOT + :PSS"+M + +e 



= PEDOT + M + : PSS~ 



(12) 



The process has referred as de-doping because the cations causes a decrease of the 
module of drain current |Ids|. This is due to the smaller number of holes available for 
conduction in PEDOT:PSS channel, as a consequence of cations incorporation into 
the PEDOT:PSS backbone. In fact, according to the Equation 12, cations adsorbed in 
PEDOT + :PSS— cause a reduction of the oxidized PEDOT+ and induces a decrease 
in conductivity upon reduction to PEDOT0. This process results reversible and when 
Vgs is switched off (Vgs = 0 V), ion diffusion occurs from the PEDOT:PSS to the 
electrolyte, increasing the number of conducting holes and, consequently, |Ids| 
(doping) 30 . 
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